Draft version December 11, 2009 

Preprint typeset using K^'T^]X style cmulatcapj v. 08/22/09 



o 
o 

(N 
o 

Q 



6 



> 

m 

O 



THERMAL TIDES IN FLUID EXTRASOLAR PLANETS 
Phil Arras^ and Aristotle Socrates^ 

Draft version December 11, 2009 

ABSTRACT 

Asynchronous rotation and orbital eccentricity lead to time-dependent irradiation of the close-in 
gas giant exoplanets - the hot Jupiters. This time-dependent surface heating gives rise to fluid 
motions which propagate throughout the planet. We investigate the ability of this "thermal tide" to 
produce a quadrupole moment which can couple to the stellar gravitational tidal force. While previous 
investigations discussed planets with solid surfaces, here we focus on entirely fluid planets in order to 
understand gas giants with small cores. The Coriolis force, thermal diffusion and self-gravity of the 
perturbations are ignored for simplicity. First, we examine the response to thermal forcing through 
analytic solutions of the fluid equations which treat the forcing frequency as a small parameter. In the 
"equilibrium tide" limit of zero frequency, fluid motion is present but does not induce a quadrupole 
moment. In the next approximation, finite frequency corrections to the equilibrium tide do lead to 
a nonzero quadrupole moment, the sign of which torques the planet away from synchronous spin. 
We then numerically solve the boundary value problem for the thermally forced, linear response of 
a planet with neutrally stratified interior and stably stratified envelope. The numerical results find 
quadrupole moments in agreement with the analytic non-resonant result at sufficiently long forcing 
period. Surprisingly, in the range of forcing periods of 1-30 days, the induced quadrupole moments 
can be far larger than the analytic result due to response of internal gravity waves which propagate in 
the radiative envelope. We discuss the relevance of our results for the spin, eccentricity and thermal 
evolution of hot Jupiters. 

Subject headings: planets 



1. INTRODUCTION 

Many of the gas giant exoplanets orbiting close to 
their parent stars - the hot Jupiters - are observed to 
have radii far larger than the radius of Jupiter, Rj, m\- 
plying high teniperat ures deep in the planetary interior 
(jBurrows et al.|[2000[ ). A powerful internal heat source 
must be present in order to prevent their rapid con- 
traction. Tidal heating is a possible solution. However, 
the gravitational tide synchronization and circularization 
timescales are much shorter than the age of the system if 
the tidal quality factor is comparable to that of Jupiter. 
Tidal heating could then power the observed radii only 
for a brief period, far shorter than the age of observed 
systems. 

We propose that time-dependent thermal forcing of 
hot Jupiters leads to a "thermal tide" torque pushing 
the planet away from synchronous rotation, and possi- 
bly circular orbits. The equilibrium spin state is set 
by the competition between the opposing thermal and 
gravitational tide torques. Gravitational tide dissipation 
continues to operate in this torque equilibrium, implying 
a steady state heat source to power the large radii. The 
thermal tide torque mechanism was initially proposed by 
Gold and Soter (1969, GS from here on) in order to ex- 
plain Venus' slow, retrograde spin. Here, we generalize 
their analysis for fluid planets in which the central solid 
core is either non-existent or dynamically unimportant. 

The existence of dynamically important thermal tide 
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torques in fluid planets has recently been cast in doubt 
(Goodman 2009 and Gu & Ogilvie 2009). These au- 
thors argue that at forcing frequencies smaller than the 
planet's dynamical frequency, the quadrupole moment 
will be negligible due to isostatic adjustment. Further- 
more, the former author asserts that the torque will push 
the planet toward synchronous rotation, rather than 
away. A careful study of the fluid motion induced by sur- 
face thermal forcing in fluid gas giant planets is needed 
to address the following questions: What is the direc- 
tion and depth dependence of the resulting flow? What 
is the relevance of the concept of "isostatic adjustment" 
to this inherently time-dependent problem? What sets 
the "dynamical frequency below which the quadrupole 
moments become small? What is the sense of the torque 
(synchronous versus asynchronous)? What is the mag- 
nitude of the torque? Once these questions have been 
answered, the ability of thermal tide torques to promote 
asynchronous rotation and tidally inflated radii can then 
be assessed. 

Here, we study the fluid flow and resulting quadrupole 
moments induced by thermal forcing at the surface of a 
gas giant planet. To isolate what we believe to be the rel- 
evant physical mechanisms, we ignore thermal diffusion, 
the Coriolis force and self-gravity of the perturbations, 
while utilizing a simplified model for the planetary struc- 
ture. 

The plan of this paper is as follows. In fj2]we provide 
motivation for performing the ensuing calculations. We 
then describe the problem set-up and details of the back- 
ground model for the fluid planet in |J21 fJHreviews the re- 
lation between density perturbation and quadrupole mo- 
ment, and the torque on the planet and orbit in terms of 
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Fig. 1. — Geometry of the semi-diurnal (m = 2) thermal and 
gravitational tidal density perturbations. The planet rotates coun- 
terclockwise. In order for the torque to push the planet away from 
synchronous spin, the density perturbation must lead the line join- 
ing planet and star. 

the quadrupole moment. The hnearized fluid equations 
are presented in 311 In ^5.11 and 15.21 we develop ana- 
lytic solutions in the zero and small, but finite, frequency 
limit. Numerical results to the linearized boundary value 
problem are presented in i j5.3l In ^ the numerical re- 
sults are explained in terms of the response of envelope 
gravity waves that are excited thermal forcing. We sum- 
marize in SjSl 

2. MOTIVATION 

We begin with a brief discussion of the possible role 
of thermal tides for the rotation rates and radii of hot 
Jupiters. 

Thermal forcing may induce a "bulge" which is mis- 
aligned with star-planet line (see figure [T]). The mis- 
alignment is a result of thermal inertia, the tendency of 
maximum temperature to lag maximum heating. Since 
insolation is finite on the day side and zero on the night 
side, thermal forcing is spread over a range of harmonics. 
In particular, the semi-diurnal (m = 2) component of the 
thermal forcing may induce a quadrupole moment in the 
fluid, which can then couple to stellar gravitational field. 
Therefore, this "thermal tidal torque" can, in principle, 
push the planet away from a synchronous spin state. The 
equilibrium spin rate is then set by a balance of the ther- 
mal tide torque, promoting asynchronous spin, and the 
gravitational tide torque, which acts to push the spin 
back toward the synchronous state. The two bulges are 
illustrated in figure [T] 

In the asynchronous spin state set by torque equilib- 
rium, the gravitational tide generates heat in a steady 
state manner. If the dissipation occurs sufficiently deep 
within the convective core of the planet, a thermal 
equilibrium will be established between (gravitational) 
tide heating of the core and heat loss at the radiative- 
convective boundary. To explain the large observed radii, 
the equilibrium core entropy must be much larger than 
that of a passively cooling planet at the same age. 

The thermal tide torque is expected to deposit angular 
momentum mainly in the surface radiative layer. Were 
there no sink of angular momentum, the total angular 
momentum of the planet would increase indefinitely. If 
we again assume that gravitational tide dissipation oc- 
curs deep in the convective core, then angular momentum 



must be transported between the thermal tide source 
region and the gravitational tide sink region. Presum- 
ably this angular momentum transport would be accom- 
plished through (turbulent) viscous forces between neigh- 
boring fluid layers. 

We estimate the importance of thermal tides in hot 
Jupiters by comparing the quadrupole moment induced 
by the thermal tide to that from dissipation acting on 
the gravitational tide. For the purposes of a fiducial 
estimate, we use the thermal tide quadrupole formula 
relevant for Venus by , Gold & Soter ( 1969, ): 
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Here AM is the mass down to the photospheric layer 
for the starlight, Rp is the planet's radius, a is the tidal 
forcing frequency and the thermal time of the absorbing 
layer is given by 
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(2) 



where Cp is the specific heat at constant pressure, T is the 
temperature, and Fi, is the stellar fiux. The frequency de- 
pendence in eq[T] reflects that the temperature and den- 
sity perturbations in the atmosphere become small when 
the forcing period is short. The accuracy with which 
eqU] represents a fluid atmosphere is the subject of the 
remainder of the paper. 

The quadrupole moment induced by the stellar gravi- 
tatio nal tide which leads to s ecular evolution is given by 
(e.g. IGoldreich fc^ote^[T96l 
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where the height of the tidal bulge is ht ^ 
Rp{Mi,/Mp)[Rp/aY , the orbital frequency is n ~ 
{GMi,/a^y^^ , the semi-major axis is a, the stellar mass 
is M+, and the tidal quality factor of the planet is Qp. 
Here a (x n — is proportional to the departure of the 
rotation rate 17 from synchronous rotation, Q = n. As a 
flducial value, the tidal quality factor of Jupiter is con- 
strained to be Qp ~ 10^ - 10*^ (Goldreich & Soter 1965; 
Yodcr & Peale 1981). 

The ratio between the thermal tide and gravitational 
tide quadrupole is given by 
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The fiducial estimate in eqd] shows that thermal 
and gravitational tidal effects are competitive for hot 
Jupiters, with large departures from synchronous spin 
a ^ n expected. 

Given the equilibrium spin rate, the gravitational tide 
will generate heat at a rate 
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This fiducial estimate of tlie heating rate is compa- 
rable, or larger t han, the core cooling rates found in 
lArras fc BildstCTi ([2006) for many "problem" planets, 
which cannot be explained by passive cooling. If this 
tidal heating is deposited sufficiently deep in the core, 
it may then prevent these planets from contracting, ex- 
plaining the inflated radii of some hot Jupiters. 

3. PROBLEM SET UP 

In this section we expand the thermal and gravitational 
forcing terms in terms of a Fourier expansion, define 
the co-ordinate system used, and discuss the background 
planetary structure model. 

3.1. description of time- dependent thermal and 
gravitational forcing 

We seek solutions to the boundary value problem de- 
scribing linear response of a gaseous planet to a grav- 
itational tidal acceleration ~'VU{x,t) and an imposed 
entropy perturbation, As(x,t), that results from time- 
dependent insolation. 

Consider a planet of mass Mp and radius Rp in a cir- 
cular orbit of separation a around a solar-type star of 
mass M^, = Mq, radius i?* = Rq and effective tem- 
perature T^, = 5800 K. The bolometric flux at the 
planet's surface is then F-^ = ashT^{R*/0')'^ , where Csb 
is the Stefan-Boltzmann constant. The orbital frequency 
is n = {GMi,/a^y^^ and the planet's spin frequency is 
n. The orbital phase in the frame corotating with the 
planet will be denoted <& = (n — il)t. 

It is convenient to express U{x,t) and As{x,t) in a 
Fourier series in time and spherical harmonics in angle 
i.e.. 



and 



As{x,t) = Asirn{r)Yi„i{ 
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The forcing frequency is then a^n = m(n — ^ For 
£ = 2, only m = ±2 are allowed due to the even parity of 
U and As. We note that while the longitude-dependent 
stellar irradiation is zero on the night side and nonzero 
on the day side, the individual Fourier components are 
nonzero on both day and night sides. By adding together 
many Fourier components one recovers the true position- 
dependent flux. 

Lagrangian and Eulerian variations are denoted by A 
and ^, respectively. The Lagrangian entropy perturba- 
tion. As, is a convenient variable to use as it varies solely 
due to non-adiabatic processes, here the time-dependent 
heating due to insolation. By contrast, the Eulerian en- 
tropy perturbations can vary due to motion parallel to a 
background entropy gradient, and the temperature can 
vary due to pressure perturbations. 

^ For an eccentric orbit, this expression is generalized to cTmfe = 
kn — where k is an integer. 



Unfortunately, the governing equations are not sepa- 
rable when the Coriolis force is included. Due to this 
complication, for simplicity, we neglect the Coriolis force. 
By doing so, the response for each individual spherical 
harmonic can studied separately. We will suppress the 
subscripts {l^ m) where no confusion will arise. 

Scalar quantities are represented as f{x,t) — 
/('')y«m(^, 0)e~*'^*, and the Lagrangian displacement 
vector, ^{x,t) is expanded in the poloidal harmonics 



^{x,t) = [^r{r) er+^h{r)rV]Yi 



(8) 



The entropy perturbation As is related to the La- 
grangian density and pressure perturbations, Ap and Ap 
respectively, via the thermodynamic relation 



Ap _ Ap 
P Tip 



As 



(9) 



Here Fi = {dlnP/d\np)s is the first adiabatic index, 
Cp — T{ds/dT)p is the specific heat per unit mass at 
constant pressure, and ps — Cp{d In p/ds)p is the negative 
of the volume expansion coefficient at constant pressure. 
For an ideal gas, ps = —1. 

If the thermal inertia of the absorbing layer is large, 
then thermal diffusion can be ignored on the timescale 
(7^^ of the perturbation. From here on, we assume that 
this limit applies. Perturbations of the heat equation 
then take on the simple form 



i Se 
a CpT 
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Here Se is the perturbed rate of heating, per unit mass. 
The depth dependence of 5e is approximated as 



6e{p) = K^,Fi, exp{-Ki,p/ g) 



(11) 



for the £ = \m\ = 2 harmonics of interest. Here, K^, is 
the effective opacity for where starlight is absorbed, at 
a characteristic pressure g/Ki,. The thermal time at the 
base of the heated layer is tth — CpT/ k-i,F-i,. Therefore, 
the imposed entropy perturbation can be written as 



— = -r- exp(-K*p/5)- 

Cp crcth 



(12) 



For semi-diurnal forcing of an asynchronously rotating 
planet, the stellar insolation Se achieves its maximum 
value at noon and at midnight. The entropy perturbation 
As then attains maxima at 3PM and 3AM, and minima 
at 9AM and 9PM. 
The £ — 2 coefficients for the tidal potential are given 

by 



U2 
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for both TO = ±2. 



3.2. Background Model 

Hot Jupiters possess a convective core, blanketed by a 
thin radiative envelope at the surf ace (e.g. iBurrows et al.l 
[2OOOI I Arras fc Bildsten|[200l) . Rather than construct a 
model with a detailed equation of state and radiative 
transfer, we use a simple parametrized equation of state 



4 



o 

^1 



10= 



10' 



1 r 



0.1 



10- 



- 1 1 1 1 1 1 1 1 1 1 1 


... 1 ... 1 . 1. - 




/ _ 






- 




.. 1 ... 1 ... 1 


... 1 ... 1 . . 



4 6 8 10 

log pressure [dyne 



cm 



12 



Fig. 2. — Propagation diagram showing Brunt- Vaisalla frequency 
{N, solid line) and Lamb frequency for £ = 2 {6^i''^c/r, dashed line) 
as a function of depth. An acoustic wave with a fixed frequency a 
can propagate in regions where a > N, ck± and gravity waves can 
propagate where a < N, ck± . Here c is the adiabatic sound speed 
and kj_ = \/E/r is the horizontal wavenumber. The region near 
the surface with A'-^ > is stably stratified, and the central region 
with A''^ ~ is neutrally stratified. 

which enforces neutral stratification in the core and sta- 
ble stratification with roughly constant scale height in 
the envelope. The equation of hydrostatic balance is in- 
tegrated using 



pip)' 



-p/pi 



1 



-p/pb 



(14) 



The constants pb, Kc and a — (phKc)^^^ are the pressure 
at the base of the radiative envelope, parameter setting 
the entropy of the core and the isothermal sound speed 
of the envelope, respectively. The factors with e^'Pl'P*' 
enforce a narrow transition at the radiative-convective 
boundary. The core compressibility is set to ~ GR^j. 
The propagation diagram for a model with p;, = 100 bar 
is shown in figure [2] The mass and radius are Alp = 
O.TMj and Rp = \.27R,j. The mass of the radiative 
envelope is roughly pb/p{r — 0) = 10^^ of the planet's 
total mass. Though this model is simple, it accurately 
reproduces the gross features of more detailed models (cf. 
Arras & Bildsten). 

To compute the sound speed we take Fi = 2 so that the 
core is neutrally stratified. The Brunt- Vaisalla frequency 
is then 
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where g{r) = Gm{r)/r'^ is the gravity, H = p/ pg is the 
pressure scale height, and the derivative ci In p/c? In p is 



computed analytically from eq. p4p . 

4. QUADRUPOLE MOMENT AND TIDAL TORQUES 

The interaction Hamiltonian, which coupled the den- 
sity field of the planet to the gravitatio nal tid al potential 
of the star is given by (e.g. i Ne wcombl 1 1 9 6 2f ) 

H= ( d^x*p*{x*,t)U{x\t). 



(16) 

Here x* and p* are the position and mass density of a 
fluid element in the perturbed planet. The fluid posi- 
tion in the background planet, x, and in the perturbed 
planet, x* , are related by a;* = a; -I- ^(a;,i), where ^ is 
the Lagrangian displacement vector. Mass conservation 
implies d^x* p*{x* ,t) = d^xp{x), where the volume inte- 
gral is performed over the background planet. Eq. (|16p 
may then be expanded as 

H = J d^xp{x)U{x + ^,t) 

d^xpix) [Uix] + I • VU{x, t) + e)(e')] . (17) 



The first term in eq. (|T7)) corresponds to the interaction 
energy of the background, and can be ignored since it 
is independent of ^. The second term is the interaction 
energy of the perturbations at linear order and is the 
term of interest. We ignore nonlinear terms of order 
on the assumption that they are small. 

To evaluate eq. (jl7p . we expand the tidal potential in 
spherical harmonics (see eqUi) 

and define the time-dependent multipolc moments of the 
planet 



Qlrnit) = J d^X p^{x, t) ■ V (/r,:„(0, (f)) ■ 



(19) 



Here D and <& are the orbital separation and phase. In- 
tegrating by parts and using the continuity equation for 
the perturbations 
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eq. (|T9|) becomes 



em{t)= / d'x r'Yl^{e,(l))Sp{x,t). 



(20) 



(21) 



Hence it is the Eulcrian density perturbation which is 
needed for the quadrupole moment. 

Since r^6p{x,t) is a real quantity, and Y^*^ ~ 
(— l)™y£_m, the moments must satisfy — 
Qe,-mi-ir- By defining We^ = [47r/(2^ + 
l)]i^m(7r/2, 0), the interaction Hamiltonian in eq. (|16p 
may be conveniently expressed in terms of a sum over 
spherical harmonics as 



(22) 



The rate of change of orbital angular momentum, 
Lorh = GMpMi,[Ga/{Mp + M^,))^/^, and the torque on 
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the planet, A/", are then 



(23) 



The secular evolution of the planet-orbit system is due 
to the presence of quadrupole moments that are out of 
phase with the tidal acceleration i.e., the imaginary com- 
ponent of the quadrupole. The torque on the planet for 
a circular orbit is dominated by the semi-diurnal term 
with |m| 2 and forcing frequency 2(n — il). The torque 
at this order is 
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(24) 



The sign in eq. ((24|l is such that density perturbations 
which lead maximum heating tend to torque the planet 
away from synchronous rotation, and vice versa. 

5. SOLUTION OF THE BOUNDARY VALUE PROBLEM 

Conservation of mass and momentum, the first law of 
thermodynamics, and two boundary conditions are em- 
ployed in order to obtain the planet's response to thermal 
and gravitational forcing. First we develop equations 
that can be solved as a boundary value problem. We 
then show that planet's response can be represented as a 
sum of cigcnmodcs. The governing equatio ns below are 
standa rd and can be found in, for example, lUnno et al.l 

Conservation of horizontal and vertical momentum are 
given by 



u 



and 



2^ 1 dSp Sp dU 
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p dr p dr 



(25) 



(26) 



respectively. Here 6p and Sp are the Eulerian pressure 
and density perturbations. We ignore the self-gravity of 
the perturbations for simplicity. The continuity eq. (j20p 
may be written as 



Sp 
P 



1 d 
r^p dr 

1 d 
r'^p dr 



(27) 



where we have used eq. ([25)) . By substituting the relation 
between Lagrangian and Eulerian quantities i.e., A = 
5 -I- 1 • V, eq. © becomes 
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(29) 



where = Tip/p is the adiabatic sound speed. The im- 
posed entropy perturbations drive vertical fluid motion 
through buoyancy forces —gppsAs/cp. The vertical fluid 
motion in turn drives horizontal motion. 

We solve the system of 3 equations ([Ml), ^ and 
for the three variables Sp, Sp and ^r- While Sp can be 
eliminated as well, we find that better numerical preci- 
sion is achieved by keeping Sp in the system of equations. 



Finally, we remark on the boundary conditions. 
Eq. (|29)) is algebraic, so no boundary condition is needed 
as this equation can be evaluated on the boundaries. At 
the center of the planet , we require all va riables to be fi- 
nite. This implies (e.g. lUnno et al.lll989l : take en r^~^ 
and Sp, Sp oc in eqO 
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(30) 



We present results for two different upper boundary con- 
ditions. The "standard" boundary condition requires 
that the Lagrangi an pressure perturbation Ap vanishes 



that the Lagrangi ; 
(|Unno et al.lll989f l 



Sp/p = 9^r 



(31) 



If the fluid perturbation is evanescent, this boundary con- 
dition is valid. We also present results for an "outgoing 
wave" boundary condition, in which wave e nergy propa- 
gates outward at the upper boundary (e.g. lUnno et al.l 
Il989f ). We implement this outgoing wave boundary con- 
dition as follows. The atmosphere above the upper 
boundary is idealized as isothermal, with constant H, N, 
g, c and horizontal wavenumber k± = \/ t(l -f \)/ Rp. In 
an isothermal atmosphere, the fluid variables Sp/p,S,r 
e^^ (cf. iGoldreich fc Kunia^ 119901 ) . where z is altitude 
and the complex constant 
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(32) 



determines the run with height. When the argument of 
the square root is positive, the wave is propagating. The 
sign of the second order wave flux is 



-^wavc — Sp ^ ^ 



r )Im(6). 



(33) 



We choose outgoing wave flux by choosing the appropri- 
ate sign of Im(s) to make F^avc > 0, given the signs of 
(7 and N'^ — a^. For the propagating case, the boundary 
condition is enforced as 



(34) 



— ( ^\ =b(— 
dr\p) \p 

When the argument of the square root in eq. ([5^ is 
negative, the wave is evanescent and we enforce the "hy- 
drostatic" boundary condition in eq. (|31[) . 



5.1. Zero- frequency limit: the equilibrium tide 

Th e equilibrium tide solution is found by setting cr — > 
(e.g. IGoldreich fc NicholsonlfT989h . The pressure pertur- 
bation is found from eq. (|25|) to be 



(35) 



jp(<=q) = -pu. 

Substitution of eq. (|35p into eq. (pS)) gives 

Sp^°°^'^ _ d\np 
P dr \g 

Inserting eg. ([55)1 and eq. ([5^ into eq. P9|) gives 

, As 

N^r^ U-gp, — . 

9 Cp 



(36) 
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Fig. 3.— 



Plot 



of the function / = 10 
K-kP/a and its first and second derivative ver- 



T 

sus depth. Here H = p/{pg) is the pressure scale heig ht. T he first 
derivative is needed for the horizontal displacement (ea l41| l. In this 
plot, the base of the heating layer is at g/re* = 1 bar and the base 
of the radiative zone is at p), = 100 bar. 

In stably stratified regions, >0 and eq. 
solved for the vertical displacement 



37t can be 



(eq)- 



U 
9 



9 



(38) 



In the central convection zone, where iV^ ~ 0, no so- 
lution exists for eq. ([37|) . since the terms with N'^ go 
to zero while the As term is nonzero. It follows that 
the concept of an equilibrium tide breaks down in re- 
gions that are neutrally stratified and therefore, the forc- 
ing frequency cannot be safely set to zero. That is, the 
fluid response in neutrally stratified regions is inherently 
a "non-equilibrium tide." Finally, substitution of eg . (I35p . 
([55]) and ([55]) into eq. (ITT)) yields the equilibrium hori- 
zontal displacement 



1 



d 
dr 



(oq) 



(eq) 



(39) 



Hence both vertical and horizontal motions exist in the 
equilibrium tide limit, driven by gravity and entropy fluc- 
tuations. 

Clearly, the equilibrium tide density perturbation in 
eq. ([55]) . and thus the equilibrium tide quadrupole mo- 
ment, has no contribution from time-dependent surface 
heating. Inspection of eq. show that the terms oc As 
cancel out one another. The last term in eq. (|29p repre- 
sents a density decrease due to heating, which is precisely 
compensated by the N'^S,r/ 9 term that results from mov- 
ing across surfaces of constant pressure. In other words, 
denser fluid is brought up from below, exactly canceling 
the local density decrease due to heating. 

Time-dependent insolation heats the surface and initi- 
ates horizontal motion. Consider subsynchronous rota- 



tion {n—Q, > 0) so that an observer at rest with respect to 
the planet sees the star rotate with angle ^{t) — (n — il)t 
in the positive (j> direction (see figure [T]). For m = 2, the 
heating function de oc cos[2(0 — $)], and maxima in the 
heating function occur at — $ — 0, tt. The entropy 
perturbation takes the form (eqfT2| 



As{x,t) = 



-1 



(40) 



n — Q J \ CpT 
xY22{e,0)sin[2{(f>~^)] 

giving the hottest points (As a maximum) at (f> — 
$ = — 7r/4, 37r/4, and the coldest points at — $ = 
7r/4, — 37r/4. This phase shift reflects the lag between 
maximum heating and maximum temperature. The dis- 
placement in the (j) direction becomes feq lSQ)) 



1 



n 



1 

3rp 



xr22(6',0)cos [2(?i-$)]. 



(41) 

The quantity in parenthesis and its derivatives are shown 
in figure m It varies mainly as as pexp{—Ki,p/g), which 
has a sharp maximum near the base of the heating layer 
at p = g/ni,. Above the base, this derivative is nega- 
tive and ^0(0) oc cos[2(0 — $)] describes motion away 
from regions of high entropy, and toward regions of low 
entropy. Were there no return flow, mass would accumu- 
late at — $ — 7r/4, — 37r/4, leading the star. For mass 
"bumps" leading the star, the quadrupole moment would 
act to decrease the planet's spin, pushing it further from 
synchronous spin. However, there is a return flow, since 
below the base the direction of ^0 reverses (see sign of 
df /dr in figure [3]). As a result, there is no mass accumu- 
lation, the density perturbation is identically zero, and 
there is no net quadrupole. What is required to drive 
asychronous spin is a net flow, integrated over depth, 
away from the the hottest points in the atmosphere. 

5.2. Finite frequency correction 

In this section we derive an analytic expression for the 
density perturbation due to thermal forcing by using per- 
turbation theory in powers of the frequency. While in- 
structive, this limit is only applicable when resonances 
with internal waves are unimportant, i.e. at forcing pe- 
riods <C 1 day or ^ 1 month. 

Let Sp and denote the complete solution including 
both the (7 = equilibrium tide and the finite frequency 
corrections. To first order in ct^, conservation of horizon- 
tal momentum gives 



dp=-pU + a^rp&'^'> = Sp'^'^i^ 



2 /-(eq) 



(42) 

Substituting eql42] into eq. ([26]) allows us to write the 
density perturbation as 



dp^Sp'^'"^^ + — 
9 



Pi. 



(eq) 



(43) 



For now, we ignore the tidal potential U and use eq. p9 
to write 



5p='- 
9 



Pi, 



(eq) 



1 



(.{t + 1) dr2 



(44) 
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This equation shows that finite fluid inertia effects, scal- 
ing as cr^, give rise to a nonzero density perturbation and 
consequently, a quadrupole moment. 

Attempts to directly integrate eq. ((44|) lead to numer- 
ical difficulties due to the second derivative term, which 
oscillates with depth. To the extent that r^^^ / g is con- 
stant, this expression gives a perfect derivative, and the 
integral would depend only on the endpoints at the center 
and surface of the planet, where the integrand is negligi- 
ble. In Appendix El we show that integration by parts 
and use of the equations of motion can be used the ob- 
tain a more monatonic integrand with significantly less 
cancellation error. 

The expression in eq. (|A3p can be estimated analyti- 
cally in the low frequency limit. If we ignore the 6p <yi 
term, substitute eq. (|38p for and treat the interior 
mass m(r) 



Mp as a constant, we find 



(4 + £)(3 + ^) 



1 



drpr 



This formula directly gives the quadrupole moment in 
terms of the applied entropy perturbation. Note that 
the ^-dependent prefactor is equal to 4 for £ = 2, and 

vanishes as £ — > OO; 

By comparison. IGold fc Soteil (|1969( ) assumed constant 
pressure and ignored vertical motion to obtain a density 
perturbation 5p/ p — —5T/T = — As/cp (for ps — — l).** 
This gives the quadrupole moment 



Q 



(GS). 



drpr + Ps 



As 



(46) 



The integrand in eq. (|45p has the same sign as that of 
eq l46[ implying a torque which generates asynchronous 
spin, but suffers a reduction relative to their expression 
by a factor 



(4-f£)(3 + £) 

e{i + 1) 



1 



(47) 



As we will see, our numerical solutions to the boundary 
value problem asymptote to the expression in eq. (|45p 
at long forcing periods. 

5.3. numerical results 

In this section we present numerical results for the so- 
lution of eq.®, ^ and ^ for Sp, Sp and ^r- To 
solve these equations, we finite difference in radius with 
second order accuracy, and write the resulting inhomogc- 
neous equation in matrix form as Mx ~ B, where M , x 
and B denote the differential operators, solution vector 
((5p, 5p and as a function of radius), and forcing vector 
(involving U and As). We find the solution of this linear 
system using a band-diagonal solver_^ The quadrupole 
moment is then evaluated using eq. (jA3p . 

Figure 0] shows the quadrupole moment Q22 (f ) as a 
function of forcing frequency a for the standard bound- 
ary condition in eq. (|3ip . Forcing is solely by As, 
i.e. U = {). The parameters used are Mp = 0.7Mj, 
Rp = l.3Rj, a — 0.05 AU, a solar-type star, base of the 
heating layer at ^/k* ~ 1 bar and base of the radiative 




It is unclear if lGold Ik Soterl II1969I ) assume (5p = or Ap = 0. 



0.1 1 10 102 

forcing period[days] 



Fig. 4. — Imaginary part of the £ = m = 2 quadrupole moment 
as a function of forcing period 27r/cr = 7r/(n — Q). The quantity 
|/m(Q22{o"))| is plotted, and short-dashed blue (solid red) lines 
show positive (negative) values. For comparison, the approxima- 
tion of Gold & Soter (1969) is s how n as the red dot-long dash line, 
while the approximation in eq. I I45I I is shown as the red dotted line. 
The peaks are due to resonances with normal modes of oscillations, 
with g-modes at long period and f-p modes at short period. The 
base of the heating layer is at pressure g/zt* = 1 bar, and the base 
of the radiative zone is at pi, = 100 bar. The solid red lines imply 
a torque pushing the planet away from synchronous spin, and visa 
versa for the short-dashed blue lines. 



zone at pb = 100 bar. The absolute value of the imagi- 
nary part is plotted. The sign is shown by the line type, 
solid red line for Im{Q22{<^)) < 0, and dashed blue for 
Im{Q22{o')) > 0. The positive forcing frequencies used 
imply sub-synchronous rotation. Torque has the same 
sign as Im{Q22), so that Im{Q22) < drives the planet 
further from the synchronous state, and vice versa. It is 
immediately apparent that the response to thermal forc- 
ing is dominated by the low order g-modes at periods 
2tt/(t ~ 1 — 30 days. This is in contrast to the gravita- 
tional tide, which has the largest response for the f-mode 
at 27r/(7 ~ 0.1 days. Consequently, the equilibrium tide 
limit does not apply until the frequency is well below the 
periods of low order g-modes. In other words, the low 
frequency limit of H5.2I is not a good approximation un- 
til periods 27r/(T > 30 days for the model shown, rather 
than 27r/(T ~ 0.1 days, as one would suspect for the grav- 
itational tide. 

For comparison, figure [4] shows the GS approximation 
in eqHB] and the finite frequency correction to the equi- 
librium tide given by eq. ((15)) . The amplitude at the 
resonances is set by the frequency spacing used to make 
the plot; for finer frequency spacing the peaks would be 
larger. The solution asymptotes to eq. (|45l) for forcing 
periods much longer than the low order g-modes. In the 
range 1-30 days, the quadrupole moment is much larger 
than the value given by eq. even midway between 

resonances. The torque alternates sign across some reso- 
nances, but not others. For short forcing periods ~ 0.1 — 1 
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Fig. 5. — Same as figure |4] but with an outgoing wave boundary 
condition at tlie upper boundary. 

day, the torque acts to synchronize the spin, reinforcing 
the dissipative gravitational tide torque. However, in be- 
tween the two lowest order g-modes, at periods ~ 1 — 2 
days, the thermal tide torque is large and acts to torque 
the planet away from a synchronous spin state. 

The sharp peaks in figure [3] are due to resonant re- 
sponse. In the absence of damping, the energy of a stand- 
ing wave diverges (in linear theory) as the resonance is 
approached. This is due to the fact that waves can re- 
flect back and forth in the resonant cavity, continually 
being pumped by the forcing. A second limit one can 
imagine is that waves with sufficiently short wavelengths 
at the top of the atmosphere (above the photosphere) 
propagate upward, carrying away their energy. Figure [5] 
shows the results for the same parameters as in figure |4l 
the only modification being the upper boundary condi- 
tion. Again, there is no forcing by the gravitational tide 
{U = 0). The outgoing wave boundary condition has 
eliminated the sharp resonances for high order g-modes. 
The exact numerical response agrees well with the an- 
alytic formula in eq. (j45p beyond a forcing period of 1 
month. For the important period range 1 day-1 month, 
the quadrupole moment is again 1-3 orders of magnitude 
larger than expected from the analytic formula, due to 
the response of low order g-modes - even for the leaky 
upper boundary. As we show in the next section, this 
enhanced response is due to the good overlap of gravity 
waves with the thermal forcing. 

6. OVERLAP OF OSCILLATION MODES WITH THERMAL 
AND GRAVITATIONAL FORCING 

At forcing periods > 1 month, the frequency depen- 
dent quadrupole moment asymptotically approaches the 
analytic result in eq. (|45|) . derived by treating fluid in- 
ertia as a small correction to the equilibrium tide solu- 
tion. However, the range of forcing periods of interest 
in determining asynchronous rotation of hot Jupiters is 
likely < 1 month. In this range of forcing periods, low 



radial order gravity waves propagate in the radiative en- 
velope and even the off-resonant response of these waves 
increases the quadrupole moment by 1-3 orders of mag- 
nitude compared to eq. (|45p . In what follows, we show 
that the forced entropy fluctuations may be understood 
as an external force that is exerted upon the planetary 
fluid. Given the expression for this force, the standard 
machinery developed for the purpose of understanding 
dynamic gravitational tides (Cowling 1941; Zahn 1970, 
1975; Press & Teukolsky 1977; Goldreich & Nicholson 
1989) can be applied to the thermal forcing case in order 
to understand the corresponding quadrupole moments. 

6.1. Forced Oscillator Formalism 

The fluid fluctuations in the planet can be thought of 
as a forced oscillator that obeys 



^ + C^= D I 



a 



(48) 



where C is a self-adjoint operator, responsible for the 
adiabatic restoring forces (Lynden-Bell & Ostriker 1967) 
and a = a*^ + is the sum of the gravitational tidal 
force 



-VU 



(49) 



and the force due to time-dependent insolation, . We 
have also added an ad hoc damping term D • ^. To com- 
pute the form of , the pressure and buoyancy forces 
in eq. (|26p must be expressed in terms of ^ using eq. [20] 
and [21] using 

C.r - Ps 

P 9 Cp 



Sp ~ pc 



PPs— ■ (50) 



V • (p^) + p — ^r 

9 

Terms involving ^ combine to give C • The total pres- 
sure force contains a non-adiabatic term that results from 
time-dependent insolation, which we identify as 



1. 



pc^p, 



As 



(51) 



The entropy fluctuation As is specified by eq. ([TT 

The externally forced fluid displacement ^ can be de- 
composed into the free adiabatic eigenmodes as 



l(x,t) 



9a(0^a(x)- 



(52) 



Here, Qait) is the time-dependent eigenmode amplitude. 
The free adiabatic eigenmodes $a{x)e~^°'''* with eigen- 
frequency ctq obey 



<«a + C-|„ = 0. 



(53) 



The Hermitian nature of the operator C implies that 
the eigenf unctions obey the orthogonality condition 
Jd'^xp^* • = Aaa^afS, whcre Sajs is the delta func- 
tion and Aaa — J d^xp\^^\'^. 

By projecting eq. (|48|l onto mode a using the orthog- 
onality relation, the amplitude qa{t) obeys the following 
forced oscillator equation 



1 



d^xpC-a{x,t), (54) 



9 



where 7q is the damping rate. In general, damping can 
be due to microphysical process such as viscosity and 
thermal diffusion. It may also represent the loss of en- 
ergy at the upper boundary due to the outgoing wave 
boundary condition, which can also then lead to a broad- 
ening of the resonant response. Plugging in a particular 
harmonic a{x,t) = a(£c, cr)e^"^*, the amplitude of the 
forced response is given by 



(55) 



Note that the mode amplitude (/^(cr) can be forced by 
both tidal gravity, which leads to an in phase response, 
and the time-dependent insolation, which leads to an out 
of phase response. 

6.2. Quadrupole Moment in Terms of Overlap Integrals 

The frequency dependent quadrupole moment is found 
by combining the mode amplitude in ea l55l and the 
quadrupole moment from eq ll9l to find 

Qlmicr) 

[J d^xpC ■ a{x, a)] [J d^xp^^ ■ V (r^F,„ 



E 



2al 



(57) 



where we identify the quadrupole moment of mode a as 
[/ d^xpC ■ a{x, a)] [J d^xp^^ ■ V {r'Y,„,{e, , 



For forcing frequencies cr <C CTq,, eg 1571 shows that 2Qaim 
can be identified with the quadrupole moment contribu- 
tion from mode a at low frequency. Note that Qaim is a 
physical quantity, since the normalization factors for the 
eigenmodes cancel out. Also note that the Qaim have 
frequency dependence solely due to As oc a~^. This de- 
pendence is factored out by writing 

Im[QQ£m(o')] =Im[(5Q£m(crQ)] ■ (59) 

The imaginary part of the quadrupole moment, which 
determines the secular torque on the orbit and spin, can 
be written out as 



lm[Qem{o-)] 



2ai 



y. 

X {Re[Qae„i{cr)]jaCra + lm[Qaim{<j)]icrl ~ (7^)) .(60) 

In the absence of dissipation (7^ — 0), the gravitational 
tide exerts no torque since the response is in phase, and 
Im[(5Q,£„i((T)] = 0. Including dissipation (70, 7^ 0), the 
gravitational tide gives rise to a torque due to the first 
term in eq. ([60]) . By contrast, the thermal tide response 
is inherently out of phase. Hence even when 7^ — 0, 
the thermal tide causes a secular torque. Inclusion of 
damping for the thermal tide will act to prevent divergent 
response at resonances, due to the denominator of the 
Lorentzian in eq. I|60p. 




0.1 1 10 

forcing period[days] 

Fig. 6. — Quadrupole moments Qa22 for the oscillation modes. 
G-modes have periods longer than 1 day, and f-p modes have peri- 
ods shorter than 1 day. The plot sho ws Q a22 for I = ra = 2 for each 
oscillation mode for the sum in eq. I I57I I for the frequency depen- 
dent quadrupole. The base of the radiative zone is at = lOObar, 
and the base of the heated layer is at g/K* = Ibar. Other param- 
eters are the same as in figure the same as in figure |4] The circles 
(squares) show modes with /m(Q„22) < {I'm{Qc,22) > 0). A ra- 
dial resolution of 9447 grid points was used to compute the results 
shown in this plot. The results are numerically converged for most 
) of the modes shown even at much lower resolution. 

Figure [S] shows the quadrupole moments Qa22i<^a) 
for each mode as defined in eq. (fSTj) and (f58|) for the 
same parameters used to find the frequency dependent 
quadrupole moment in figure [H These quadrupole mo- 
ments were computed by fitting the (522(c) near reso- 
nance using eq. (|57p in order to read off the coefficient 
Qa22 (fct). It is immediately apparent that the domi- 
nant overlap is with the lowest order g-mode at period 
~ 1 day. The f and p mode quadrupoles are down by 1-2 
orders of magnitude. 

To verify the equivalence between the solution of the 
boundary value problem and the sum over eigenmodes, 
the frequency dependent response (322(0") is compared 
to the sum over eigenmodes from eq. ([ST)) in figure [71 
Damping is neglected, i.e. the standard boundary con- 
dition in eq. (PT|) is used. Despite the fact that only a 
handful of low order g-modes are used in the sum, and 
that the f and p modes are not included, there is good 
agreement with the solution to the full boundary value 
problem. 

The difference resulting from the two different up- 
per boundary conditions as shown in figures H] and [5] 
can be understood using eql60l Figure [4] has 7^ = 0, 
and hence divergences appear in the quadrupole mo- 
ment near the frequencies of internal standing waves. 
In figure EJ damping is introduced through the out- 
going wave boundary condition. In fact this damping 
is so large that 7^ ~ CTq, and the Lorentzian factors 
o'Q/((o'a -o-^)^+7aCr2) =0(1). For Lorcuztian factor of 
order unity, the induced quadrupole moment is of order 
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Fig. 7. — Comparison of the nu meri cal r esp onse (black lines) 
versus the sum over modes in eq. 1571 and 1591 (red lines) . The 
coefficients Qa22{'^a) for each mode are shown in figure [6l and the 
response is for the case shown in figure|4] Only modes in the period 
range 1-5 days were used to make the red lines; in particular, no f 
or p modes were used. Note the sum over modes agrees well for low 
order g-modes, even midway between resonances, implying the off- 
resonant response of the f and p-modes is weak compared to that 
of low order g-modes. Agreement could be improved by including 
additional modes in the sum. 



the overlap integral Qaimio'a)- 
Even for large damping rates 



a a, the low order 



g-modes still dominate the response since their overlap 
integrals with the thermal tide force are far larger than 
the f and p-modes (see figure [6|) . In the duration of one 
wave travel time in the radiative layer, the wave can build 
up enough energy to dominate the response. The overlap 
integrals in figure [S] are larger than near-equilibrium tide 
result from eq. (j45[) in the period range 1 day - 1 month, 
precisely the region where enhanced response is seen in 
figure m 

The spatial profiles of thermal and gravitational forc- 
ing differ greatly from one another. Thermal forcing 
is localized in physical space to the outer layers of the 
planet, resembling a delta function in space. By contrast, 
the gravitational tide most closely resembles the f-mode, 
and ca n be thought of as a delta function in momentum 
space eiseneggeilll994[ ). In other words, since thermal 
forcing is narrow in physical space, the fluid response, 
in terms of normal modes, is broad in momentum space, 
and visa versa for the gravitational tide. However, the 
modes which propagate within the stably stratified en- 
velope dominate since their frequencies are comparable 
to forcing frequencies of interest, and their eigenfunctions 
are confined to the radiative layer where the thermal forc- 
ing occurs. The frequencies of f and p-modes are too high 
to be resonant, and their eigenfunctions extend over the 
entire planet. 

7. ON THE NEGLECT OF THE CORIOLIS FORCE, 
THERMAL DIFFUSION, AND BACKGROUND FLOWS 



In an attempt to elucidate the origin of quadrupole 
moments from thermal forcing, we have made a number 
of simplifying assumptions so as not to obscure the basic 
physics. In this section we comment on how our results 
might be changed, both qualitatively and quantitatively, 
by relaxing these assumptions. 

7.1. the Coriolis force 

Consider the effect of uniform rotation, with spin fre- 
quency fl. When a < ft, the Coriolis force becomes 
important for the fluid response and we anticipate sev- 
eral additional effects due its presence. First, the gravity 
wave dispersion relation and fluid motions are altered. 
For a ~ f2, we expect the change in the quadrupole to 
be at the factor of a few level in the "bump" between 1-30 
days. Second, new wave families arise, mainly restored 
by the Coriolis force - the Rossby waves, or inertial 
waves. These waves can propagate throughout both the 
radiative envelope and the convective core. We expect 
that, for sufficiently rapid supersynchronous rotation, 
Rossby resonances may give rise to a larger quadrupole 
moment. In the perturbation theory calculation of ^5.2\ 
we expect larger quadrupole moment when the spin is 
nearly synchronous, since the Coriolis terms oc afl will 
be larger than the inertia terms oc cr^. In figures S] and [5] 
this will lead to the "analytic" line becoming flat, instead 
of decreasing to longer forcing period. 

For sufficiently large quadrupole moment, the equilib- 
rium spin will be asychronous enough that a > fi, and 
we expect our results to be approximately quantitatively 
correct. Substitution of the quadrupole moment in figure 
[5] into eq. SI for Qp — 10^ we find cr ~ n in the period 
range of ~ a few days, implying the spin frequency, or- 
bital frequency and forcing frequency are all comparable 
in torque equilibrium. In this event, the Coriolis force 
is important, but will only change the frequencies and 
eigenfunctions of the gravity waves at the order unity 
level. 

7.2. thermal di 



The heating function in ea llll describes the law of ex- 
ponential attenuation of stellar radiation. By setting this 
to occur at a pressure ~ 1 bar, we are assuming there 
are no absorbers high in the atmosphere. Until recently, 
this assumption was standard, but there may now be 
evidence of an inversion layer due to such absorbers i n 
a fraction of observed hot Jupiters (|Fortnev et al.|[2008l ). 
In that case, incorpor ating thermal forcing by the bound- 
ary condition used in lGu &: Qgilvid (|2009f ) may be more 
appropriate. ^ 

The entropy perturbation in eq llOl assumes that the 
forcing period is much shorter than the thermal time in 
the heated layer, so that radiation diffusion may be ig- 
nored. In fact, thermal diffusion and damping of down- 
w ard propagatin g waves is expected, as in the solutions 
of iGu fc Qgilvid ((20091 ). Furthermore, a shorter diffu- 
sion time may partially "erase" temperature perturba- 
tions, reducing the buoyancy force upon which gravity 
waves rely. This may cause wave damping or reflection 
near the optical or infrared photospheres. By choosing 

^ If one considers a two-frequency, two-angle approximation 
to the radiative transfer, both the "greenhouse" and the "strato- 
sphere" cases may be included. 
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the "perfect reflector" and "leaky" upper boundary con- 
ditions in section [5?3l we attempted to simulate two ex- 
treme limits of the effects of radiative diffusion. While 
thermal diffusion may decrease the amplitude of gravity 
waves near the photosphere, deeper in the envelope the 
thermal time becomes sufficiently long so that damping 
can be ignored. 

For forcing period of order a few days, the thermal 
time at the base of the heating layer is of order the forc- 
ing period. Therefore, the approximation that the layer 
possesses a large thermal inertia is only approximately 
valid. Another effect we have ignored is that nonlin- 
ear effects become important when the forcing period is 
longer than the thermal time. 

7.3. fluid flow in the background state 

We now comment on the possible role of fluid motion 
such as zonal flows in the background state. Fluid motion 
in the background state does not preclude the existence 
of the gravity (or inertial) wave response as described 
in this paper. Ignoring shearing of the velocity field, 
such a flow would simulate uniform rotation, so that 
the Coriolis force becomes important. However, when 
there is velocity shear in the background state, we ex- 
pect the wave motions to be significantly modified when 
the shear rate is comparable to the wave frequency. Since 
inertial wave response relies on angular momentum gra- 
dients, we expect these waves to persist in the presence 
of velocity shear. Most atmo spheric circulation simula- 
tions for the hot Jupi ters (e.g. IDobbs-Dixon fc Liiill2008l 
iShowman et al.l[2009l ) assume synchronous rotation, and 
find powerful advective flows. Relatively few simula- 
tions with asychronous rot ation have been performed, 
but the recently results by IShowman et al.l (120091 ) are 
useful in this context. They found that when the degree 
of asynchronous spin was increased, the flow velocities 
were decreased. It remains to be seen what winds are in- 
duced for rotation rates appropriate for the torque equi- 
librium described in section [21 All current studies of at- 
mospheric circulation in hot Jupiters ignore the thermal 
tide torques. From the analysis presented in this work, 
it seems likely that the qualitative and quantitative out- 
come of these simulations will change once thermal tidal 
torques are included. 

8. SUMMARY 

We considered the ability of a planet, subjected 
to time-dependent stellar irradiation, to develop net 
quadrupole moments. The existence of such quadrupole 
moments allows the stellar tidal field to exert torques on 
the planet, possibly creating asynchronous spin. Such 
asynchronous spin could lead to large tidal heating rates, 
through the dissipative gravitational tide, perhaps suffi- 
cient to power the large observed radii. 

In tj5.ll and 15. 2[ the fluid equations are solved treating 
inertia as a small parameter. At zeroth order in the forc- 
ing frequency, a circulation pattern is found from hot 
to cold at small depths p < g/K^,, and a return flow 
from cold to hot at larger depth p > g/^*- There is no 
density perturbation associated with this flow and there- 



fore, zero quadrupole moment. The inclusion of finite 
inertia as a small perturbation gives rise to a finite den- 
sity perturbation and therefore, quadrupole moment. In 
the assumed limit of small fluid inertia, the phase of the 
quadrupole is found to have the correct sign to induce 
as ynchronous ro t ation in the planet, similar to the work 
of' lGold fc Soted (|l969l) . However, the magnitude of this 
quadrupole is reduced by a factor ~ A{a/N)'^. A propa- 
gating wave-like response of any form is eliminated in this 
"near-equilibrium tide" approximation. The value of the 
resulting quadrupole moment is insensitive to boundary 
conditions and depends solely on the local forcing in the 
atmosphere. Numerical calculations presented in figure 
m and [5] confirm that the quadrupole moment approaches 
this analytic limit at long forcing periods. 

In t j5.3i the fluid equations are solved as a boundary 
value problem. Two different upper boundary conditions 
are used, which effectively treat the radiative envelope 
as a perfect resonant cavity, or as a cavity with an open 
upper lid, allowing waves to be radiated upward. Us- 
ing both boundary conditions, which may bracket the 
true result including radiative transfer effects, it is found 
that the quadrupole moments in the period range 1 day 
- 1 month are far larger than the non-resonant near- 
equilibrium tide calculation in eq. [45] by 1-3 orders of 
magnitude. The reason for this enhanced response is 
discussed in ij6l Unlike the gravitational tidal forcing, 
time-dependent thermal forcing and low order envelope 
g-modes couple well due a favorable spatial overlap. 

In section [21 we argued that the mag nitude of the 
thermal tide torques using the formula of iGold fc Soted 
()1969( ) is sufficient to generate large deviations from syn- 
chronous rotation as well as large tidal heating rates, 
easily sufficient to power the observed planet radii. Our 
numerical results in §5.3l show that the full numerical so- 
lutions for the torque can indeed have the correct sign 
to generate asynchronous spin for some ranges of forcing 
(rotation) freque ncy. The magnitud e of the torques is 
smaller than the 'G old fc Soterl (|1969f ) result by a factor 
of a few, but still large enough to be interesting. We ex- 
pect that the inclusion of the Coriolis force will increase 
the torque at long forcing periods, since the Coriolis force 
is larger than fluid inertia in that limit. 

In future investigations, we plan on including the ef- 
fects of the Coriolis force and thermal diffusion in the 
calculation of the quadrupole moments. We expect that 
both physical effects could lead to qualitative and quanti- 
tative changes to the non-resonant and resonant thermal 
tidal response. The resulting torques will then be used 
to solve for the equilibrium values of the planet's spin, 
radius and eccentricity. 
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APPENDIX 

ACCURATE EVALUATION OF THE QUADRUPOLE MOMENT: INTEGRATION BY PARTS 

As discussed in section [521 direct integration of r^dp to find the quadrupole moment is subject to large cancellation 
errors. This behavior is somewhat surprising, and does not occur for the gravitational tide quadrupole moment. 
The origin of the problem is that 5p is proportional to the second derivative of a certain quantity, resulting in large 
cancellations during the integration. We noticed this cancellation error during numerical resolution studies in which 
the results did not converge rapidly with increasing resolution. 

Here, we show that this cancellation error can be ameliorated through successive integration by parts, and use of 
the equations of motion. By using the exact equations of motion, this method does not introduce approximations. 
Integration by parts leads to an integrand that is much more monatonic, leading to better convergence. 

The quadrupole moment is given by 



By substituting for 5p from eq. 



261) gives 

Q = 



Q 



9 



2^ 1 dSp dU 
p dr dr 



(Al) 



(A2) 



Integrating the pressure gradient term by parts, discarding the (small) boundary terms, and solving eq l28l for 6p, we 
find 



Q 



g \ dr 



dr— 

dr 



r^Sp 



£{£ + !) 



Integrating the d/dr{r'^ p£^r) term by parts and simplifying gives the final result 



Q= drr^+'^Sp^'^'i^ 



dr 



£{£+!) dr^ 



„2+e 



£{i+l) dr 



„2+l 



pU 



Sp 



(A3) 



The first term in eq. (jA3[) is just the equilibrium — quadrupole from gravitational forcing alone - see eq. p6p . 
Although the Sp term in eq. (jA3[) may be subject to cancellation error, this term should be smaller by a factor a^r/g 
compared to the other terms after integration. Since the equilibrium tide component of is monatonic, the terms 
involving do not suffer cancellation error. We find numerically that integrand is far more monotonic than r^~^^5p, 
and is less subject to cancellation error. Throughout, we use eq. (jA3p when numerically evaluating the quadrupole 
moment that results from thermal forcing. 
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